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Abstract 



The transition in random fiber networks from two-dimensional to three-dimensional 
planar structure driven by increasing coverage (total fiber length per unit area) is 
studied with a deposition model. At low coverage the network geometry depends on 
the scale-free product of fiber length and coverage while at high coverage it depends on 
a scale-free combination of flexibility, width and thickness of the fibers. With increasing 
coverage the roughness of the free surface decouples from the substrate, faster when 
fibers are stiffer. In the high coverage region roughness decreases exponentially with 
increasing fiber flexibility. 



PACS numbers: 81.30Dx, 81.15Lm, 64.60-i 
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The mechanical and structural properties of two-dimensional random systems have 
recently been investigated intensively using computational and theoretical methods 
[l], [J. However, experimental studies are still few in numbers as suitable model systems 
seem to be lacking. Planar materials with random fibrous structure, such as polymer 
films, short-fiber composites, paper and nonwovens, provide promising possibilities for 
experimental work ||, in addition to their intrinsic technological interest. Fibers in real 
materials, such as paper, are much longer than the sheet is thick so that the in-plane 
alignment of fibers is over 90% ||]. The structure is therefore planar, or "layered" (but 
without well-defined layers), and the mechanical properties are essentially those of a 
two-dimensional system. The statistical geometry of strictly two-dimensional random 
fiber networks is known quite well ||, but not the geometry of planar fiber networks. 
Analytical models |6|, [7] have so far not dealt properly with the bonding of fibers that 
are not located in adjacent "layers". This bonding is controlled by the random pore 
structure of the intervening " layers" . 

In this Letter we employ a simple computer simulation model to study how the geom- 
etry of a random network of flexible fibers forms during the build-up or sedimentation 
on a two-dimensional substrate ("growth in 2+1 D"). If coverage (i.e. the total fiber 
length per unit area) is low, the fiber networks is almost two-dimensional and the in- 
plane mechanical properties are governed by coverage. At high coverage, instead, the 
in-plane properties depend on the three-dimensional "bulk" density. However, in both 
cases the intrinsic property is the average number of connections (inter-fiber bonds) 
per unit fiber length or "coordination number" ||. This in turn depends on the dimen- 
sions and flexibility of fibers. We concentrate on the coverage- driven transition from 
2D to 3D network geometry with increasing coverage and on the corresponding change 
in the surface roughness 0. Our results show the connection between the strictly 
two-dimensional random fiber networks and their planar three-dimensional real-world 
counterparts, and thereby allow the comparison of future experiments with theoretical 



or numerical studies [10], O, 12 



We start with a two-dimensional square lattice of linear size L — 10 . . . 1000. The size L 
affects properties only close to the percolation threshold and, for surface roughness, at 
very high coverages. Fiber length If is an odd integer multiple of the lattice constant 
(= 1), fiber width and thickness are fixed, Wf = tf = 1, and emphasis is given to 
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If 3> 1 characteristic of fibers Jf3|- Fibers are aligned at random [|14[] in the two 
principal directions and the centers of mass are located, again at random, over cell 
centers. Hence local coverage c (the total fiber length per cell) has an integer value. 
Bending flexibility Tf of fibers is defined by limiting the change in the vertical elevation 
of a fiber from one lattice cell to the next to < Tf (Fig. 1). A high Tf gives a dense 
network, a low Tf a sparse one. The so-called wet fiber flexibility [|15| (WFF) of paper- 
making fibers is related to Tf through Tf = [C x Wf x WFF] 1 / 4 where the constant C 
depends on experimental details. 

During the buildup, fibers are placed down independently, one after another, as if they 
were sedimenting from a dilute suspension. Each fiber is kept straight and parallel to 
the substrate until it makes the first contact with the underlying network. Thereafter 
the fiber is deformed so as to lie as low as possible - but not below the substrate - while 
still obeying the deflection constraint Tf. Periodic boundary conditions are applied 
except in the percolation study. The discretization of network structure should not 
affect generic properties. The cumulative build-up of the network is a more serious 
simplification as it does not distinguish between densification mechanisms such as hy- 
drodynamic drag forces and pressing of a filtrated fiber mat. However, the finite width 
of real fibers implies that they tend to follow a staircase pattern like that given by the 
model (cf. Fig. 1). The construction also leads to two-sided surface properties, but 
here only the free top surface is studied. 

The model networks are characterized as follows. Local thickness t is the distance 
from the free top surface to substrate. Thus the mean thickness i includes the voids 
between the network and substrate. However, in calculating density p these voids are 
subtracted away and so p < 1 counts only for the voids within the sheet. This can 
be considered as the "true" density, one that is independent of surface roughness. A 
"coarse-grained" or "apparent" density is also employed to mimic the finite in-plane 
resolution of real measurements. In this case the system is first divided into subgroups 
of 10 by 10 cells (the size being arbitrary) and the largest value t = t g is determined 
for each subgroup. The "apparent" density is then given by p g = c/i g where c and t g 
are the mean values of c and t g . Finally, the coordination number n is defined as the 
average bonded surface area of a fiber, divided by 2lfWf |§. 
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At low mean coverage c the mechanical properties of a planar network should be con- 
trolled by c— c c where c c is the percolation threshold fT6| , |P7| for mechanical connectivity. 
We determined c c numerically using the following rules: 1) the network is connected 
when there exists a cluster of fibers that is connected to at least one cell on two oppo- 
site boundaries of the system; 2) two cells are connected by a fiber if the fiber covers 
the both sites; and 3) two fibers covering a cell are connected, if the top surface of the 
lower fiber is in contact with the bottom surface of the upper fiber. The threshold c c is 
then equivalent to the bond percolation density of the dual lattice for bonds of length 
If — 1. Thus an overall c c ~ 1/(7/ — 1) behavior is expected. On the other hand, if the 
network had continuous spatial symmetry and the fibers had zero width and infinite 
flexibility, then the percolation density would be given by c c = 5.71/// [[18]. 

The computed percolation threshold c c is shown in Fig. 2 for the lattice size L = 200. 
With long fibers (If > 20) the threshold is c c ~ 8/(7/ — 1) when Tf — > oo, which is 
slightly larger than with continuous symmetry. With short fibers the non-zero width Wf 
lowers the asymptotic value so that c c ~ 6/(7/ — 1) when // = 5. The asymptotic values 
at Tf — > oo are independent of L, increasing L upto 800 revealed no deviation from the 
behavior shown in Fig. 2. A random walk argument suggests that in the opposite limit 
of Tf — > the radius r of a cluster with N fibers is r ~ lf\/N and therefore the thickness 
of a percolating cluster is ~ tf(L/lf) 2 . Thus the network becomes three-dimensional 
and c c very large and dependent on If and L. The non-smooth small-T/ behavior of 
c c seen in Fig. 2 when If = 5 is reproducible and hence probably an artefact related 
to the fact that connectivity benefits from rational values of Tf. Since the maximum 
vertical deflection per unit length of a fiber is Tf, crossover between the flexible fiber 
and stiff fiber behaviors should occur at Tf ~ tf/lf \ our data suggests that this occurs 
at Tf « A/If. 

When coverage c increases, coordination number n increases too and eventually be- 
comes a constant, see Fig. 3. In the high-coverage region n is independent of If. In 
contrast to n, density p decreases with increasing coverage because at c — > there are 
no intra-sheet voids. An order of magnitude higher coverage than with n is needed 
before "bulk" values of p, independent of c, are obtained. The difference is especially 
clear with stiff fibers. As n should be the main factor controlling the in-plane mechan- 
ical properties of the network, it follows that these become homogeneous, independent 
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of coverage, at much lower coverages than density does. 



The length of a fiber If affects the above results at low coverages, c ~ c c . For example, 
at the percolation threshold and with large Tf the number of inter-fiber bonds per 
fiber, given by 2nlf is 7.1,6.3 and 4.7 for If = 21,11 and 5. However, at higher 
coverages the three-dimensional network structure should be controlled by the scale- 
free combination WfTf/tf. This can be rationalized by observing that each fiber on 
the surface of the network prevents new fibers from making contact with others on 
(wf + lf)tf/Tf « Iftf/Tf sites (on the average) when If 3> Wf, while bonding with the 
fiber in question is possible in Wflf cells. Measured values of the wet fiber flexibility 



T5H and known dimensions of paper-making fibers yield WfTf/tf = 0.5 — 0.7. 



Coordination number n is tedious to measure directly and therefore density is some- 
times used to estimate it. Indeed, Fig. 4 shows that if c is held constant and Tf is 
varied, then n is roughly proportional to p and p g , as long as p < 0.6 and c is not too 
low. For example, in ordinary paper p = 0.2 — 0.5 so that the linear behavior should 
be valid. The reproducible discontinuities at n « 0.8 and 0.9 (corresponding to Tf = 2 
and 4) are probably again artefacts of the discrete model. It is interesting to note that 
the slope dn/dp ~ 1.2 with > 10 and approaches unity only for short fibers. The 
slope depends slightly on coverage and thus density usually cannot be used to compare 
coordination numbers if coverage varies. 

The two-dimensional projection of the network obeys Poisson distribution || and there- 
fore the standard deviation of local coverage should be a c = y/8. Our simulations follow 
this relationship very accurately. The flat substrate affects the network structure at 
small c but at sufficiently high c the structure is truly three-dimensional, independent 
of c. Then, for example, the standard deviation of the number of pores per cell p 
becomes a p = P(T/)v / B. The amplitude PiTf) decreases with increasing Tf. Areas of 
high coverage have few pores and therefore P < 1. 

Simulation results for the roughness at of the free surface (i.e. the standard deviation 
of t) can be compactly expressed as 



2 



cx[A(T f ) + B(c)} (1) 



(for general tf and Wf we would of course expect that a t ~ tf and A = A{wfTf /tf)). 
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The functions A and B are shown in Fig. 5 for system size L = 100 and c < 10 4 . They 
are independent of fiber length as long as If > 5. Equation (1) describes the coverage- 
driven cross-over from a 2D to 3D network structure. Hence, as coverage increases 
roughness is no longer affected by the flat substrate. This corresponds to B becoming 
smaller than A. The decoupling is fast with stiff fibers (A is large) because there are 
many pores in the network. The simulated roughness values can easily be decomposed 
according to Eq. (1) since a t becomes constant (or A — > 0) at large Tf. 

The form of Eq. (1) arose from the observation that at is controlled by high "peaks" 
and deep "valleys" of the surface. Hence the variance of t could be a sum of two terms 
especially since the "peaks" are dominant at low coverage or high fiber flexibility and 
the "valleys" at high coverage or low flexibility. The exponential decay of roughness 
at a fixed c (cf. A(Tf) in Fig. 5) then follows from the decrease in the depth of 
the "valleys". The planar fiber network differs from many growth models || in that 
there is nothing that would play the role of surface tension, although fluctuations of 
surface height are limited by Tf over the fiber length. Clearly, making Tf small does 
not give a smooth surface. In fact our model is similar to the Void model |19|, |9| . The 
growth behavior at high coverage, including the effects of system size and the naturally 
included overhangs, is beyond the scope of this study and will be presented elsewhere. 

Summarizing, we have presented a computer simulation study of planar random fiber 
networks. The results can be compared with experiments on disordered fibrous ma- 
terials such as paper. At high coverage c the network is three-dimensional and its 
geometry is controlled by the scale-free combination of fiber width, thickness and flex- 
ibility, (wf/tf)Tf, provided that Wf/lf <C 1 which usually holds fibers. At low c the 
network is two-dimensional and its properties depend on the scale-free product of cov- 
erage and fiber length, c If, provided that (lf/tf)Tf is large (> 4). Thus fiber length is 
important only at low coverages. Our results connect the properties of realistic planar 
fiber networks with those of ideal two-dimensional networks. In particular, if coverage 
is held constant then coordination number is linearly related to density for p < 0.6 but 
the slope depends slightly on c. Roughness of the free surface decreases exponentially 
with increasing Tf. Because of its simplicity our model could also be used to simulate 
other properties such as three-dimensional permeability and transverse compressibility. 



7 



This study was financially supported by the Technology Development Centre of Finland 
and by the Academy of Finland which is gratefully acknowledged. 

References 

[f] Permanent address. 

[1] Statistical Models for the Fracture of Disordered Media, H.J. Herrmann and 
S. Roux, eds. (North- Holland, Amsterdam, 1990). 

[2] Sahimi, M., Physica A 186 (1992) 160. 

[3] Kertesz, J., Physica A 191 (1992) 208. 

[4] Radvan, B., in The fundamental properties of paper related to its uses, F. Bolam, 
ed. (British Paper and Board Ind. Fed., London, 1976) pp 137-147. 

[5] Corte, H. and Kallmes, O.J., Tappi J. 43 (1960) 737. 

[6] Kallmes, O.J., Corte, H. and Bernier, C, Tappi J. 44 (1961) 519. 

[7] Gorres, J., Sinclair, C.S. and Tallentire, A., Paperi ja Puu 71 (1989) 54; Gorres, J. 
and Luner, P., J. Pulp Paper Sci. 18 (1992) J127. 

[8] For paper the quantity corresponding to coordination number is relative bonded 
area, RBA, of the fibers. See Kallmes, O.J., Corte, H. and Bernier, C, Tappi J. 
46 (1963) 502. 

[9] Krug, J. and Spohn, H., in Solids far from Equilibrium: Growth, Morphology, 
Defects, C. Godreche (Cambridge University Press, Cambridge, 1991) 479 - 582. 

[10] Alava, M.J. and Ritala, R.K., J. Phys.: Condens. Matter 2 (1990) 6093; 
Alava, M.J. and Ritala, R.K., Physica Sripta T33 (1990) 155. 

[11] Astrom, J.A. and Niskanen, K.J., Europhys. Lett. 21 (1993) 557. 

[12] Astrom, J., Saarinen, S., Niskanen, K. and Kurkijarvi, J., J. Appl. Phys. 75 (1994) 
2383. 

8 



[13] Fiber thickness is arbitrary in our study. For paper-making fibers lf/tf > lf/wf ~ 
20 - 50. 

[14] We use the R250 random number generator, Kirkpatrick, S. and Stoll E., J. Comp. 
Phys. 40 (1981) 517. 

[15] Steadman, R. and Luner, P., Papermaking raw materials, V. Punton, ed. (Me- 
chanical Engineering Publ. Ltd., London, 1985) 311-338. 

[16] Kirkpatrick, S. Rev. Mod. Phys. 45 (1973) 574. 

[17] Schwartz, L.M., Feng, S., Thorpe, M.F. and Sen, P.N., Phys. Rev. 32 (1985) 4607. 
[18] Pike, G.E. and Seager, C.H., Phys. Rev. B 10 (1974) 1421. 
[19] Void, M.J., J. Colloid. Sci. 14 (1959) 168. 



9 



Figures 



FIG. 1. Two crossing fibers on the square lattice. The height of the "steps" 
is = T f = 1/3. 

FIG. 2. The percolation threshold c c versus Tf for fiber lengths If = 5, 11 
and 21 (diamonds, crosses and squares, respectively); system size L = 200. 

FIG. 3. Coordination number n, density p and "apparent" density p g (tri- 
angles, squares and diamonds, respectively) against coverage c on log-scale 
with Tf = 1 and 2 (filled and open symbols, respectively). Arrows indicate 
the percolation point. System size L = 1000, fiber length If = 21. 

FIG. 4. Coordination number n against density p and p g (open and closed 
symbols, respectively) with growing fiber flexibility Tf. Coverage c = 10.5, 
31.5 and 94.5 (triangles, squares and diamonds, respectively). Fiber length 
If = 21, system size L = 1000. 

FIG. 5. Function A{Tf) on semi-log scale (Eq. (1)), for coverages c = 
10, 100, 10 3 and 10 4 (squares, crosses, triangles and diamonds, respectively). 
Fiber length If = 5, system size L = 100. Inset shows function B(c) on 
semi-log scale. 
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